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BACKGROUND 


1 . 


Finite  element  structural  analysis  is  gaining  wide  acceptance  in  the 
field  of  projectile  design.  In  the  past,  computer  memory  and  storage 
requirements  have  limited  the  application  of  finite  element  analysis  to 
large  main  frame  computers  and  workstations.  Recently,  great 
improvements  have  been  made  in  the  computing  power  of  desk  top  personal 
computers.  As  a  result,  interest  has  developed  in  applying  finite  element 
software  to  the  P.C  environment.  The  mathematics  and  computer 
algorithms  for  finite  element  analysis,  in  particular  the  extensive  matrix 
operations  required  to  solve  the  stiffness  and  deformation  of  the 
structure,  have  existed  for  many  years,  and  are  readily  obtained  from 
engineering  references.  However,  the  most  tedious  and  time  consuming 
operation  is  the  generation  of  the  finite  element  model.  A  suitable  mesh 
generator  is,  therefore,  required  to  fully  implement  any  finite  element 
package  on  a  personal  computer.  This  report  presents  an  approach  to 
developing  the  two-dimension  finite  element  grid.  The  computer  code 
listing  is  provided  in  the  appendix,  along  with  an  example  mesh 
generation.  The  documentation  is  sufficiently  detailed  to  permit  the 
programmer  to  tailor  the  mesh  generation  software  to  suit  particular 
finite  element  analysis  requirements. 
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2 .  MESH  GENERATOR  REQUIREMENTS 


Advanced  projectile  structural  design  typically  involves  the 
analysis  of  complicated  projectile  shapes,  such  as  the  sabot,  shown  in 
cross-section,  in  Figure  2.1  Figure  2.2  shows  a  generic  long  rod 
penetrator  within  this  structural  sabot.  This  type  of  projectile 
configuration  is  typically  used  in  kinetic  energy  anti-tank  projectiles. 
The  cross-section  drawing  shown  in  Figure  2.2  is  created  using  a  generic 
solids  modeling  algorithm.  Table  2.1  shows  the  structural  dimensions 
defining  this  generic  projectile. 


Table  2.1 

Projectile  Geometric  inputs 


* 


SABOT 

2.0 

a 

1.0000 

4 

ELEN 

01 

02 

L 

Rho 

FORM 

R 

H 

K 

1 

1.3 

2.0 

2.25 

3 

2 

2.0 

2.5 

0.25 

43 

0.5 

0.2057 

1.4557 

3 

3.0 

3.0 

1.0 

2 

4 

2.5 

2.0 

0.25 

44 

0.5 

0.2057 

1.4557 

5 

2.0 

1.4 

2.25 

3 

6 

1.4 

1.65 

0.5 

43 

1.0 

0.0157 

1.70 

7 

1.6S 

3.0 

1.0 

3 

8 

3.0 

3.0 

0.5 

2 

9 

1.0 

1.0 

6.5 

2 

10 

1.0 

1.5 

0.25 

42 

.5 

0.2057 

-0.2943 

11 

1.5 

2.7 

0.85 

3 

12 

* 

2.7 

2.7 

0.40 

2 

PEN 

0.0 

1 

0 

6.6400 

ELEN 

01 

02 

L 

Rho 

FORM 

R 

H 

K 

1 

1.0 

1.0 

10.5 

2 

** 
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Figure  2.1 

Cross-Section  of  a  Generic  Profectiie  Sabot 


Figure  2.2 

Cross-Section  of  Generic  Kinetic 


Energy  Projectiie 
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What  is  desired  is  to  create  a  finite  element  grid  within  the 
boundaries  of  this  projectile  form.  Finite  element  analysis  typically 
requires  the  structure  to  be  divided  into  individual  quadrilateral  elements. 
By  convention,  each  element  is  defined  by  four  corner  nodes,  in  either  a 
clock-wise  or  counter  clock-wise  direction,  depending  on  the  finite 
element  algorithm.  To  ensure  that  all  of  the  elements  are  structurally 
connected,  adjacent  elements  are  defined  by  common  node  numbers, 
depending  on  their  location  in  the  assembled  grid.  A  file  structure  is  then 
developed  which  allows  the  finite  element  analysis  program  to  create  the 
stiffness  matrix,  based  on  the  element  corner  nodes  and  the  x  and  y  or  z 
and  r  coordinates  of  each  node. 

3.  CREATING  THE  FINITE  ELEMENT  GRID 

The  mesh  generator  algorithm  presented  in  this  report  follows  a 
systematic  process  to  creating  the  finite  element  grid.  First  the 
projectile  form  is  subdivided  into  homogenous  areas,  by  defining  area 
corner  points  and  lines,  based  on  the  modeler’s  vision  of  how  the  grid  is  to 
be  developed.  Figure  3.0  shows  an  outline  of  the  structure,  when  broken 
into  sub-areas  for  grid  generation.  Figure  3.1  shows  32  discrete  points 
defining  the  boundary  of  this  structure.  Table  3.1  shows  the  input  format 
for  these  32  points,  as  required  by  the  software.  Four  additional  points 
are  entered  in  this  example,  which  are  not  corner  points,  but  are  required 
for  complete  description  of  the  geometry.  The  structure  in  Figure  2.2 
includes  four  radii,  which  describe  the  curvature  of  the  outline  at  the  base 
of  the  central  bulkhead  and  the  root  of  the  front  scoop.  For  each  radius  in 
the  structure,  the  center  point  of  the  curvature  also  needs  to  be  defined. 
Therefore,  points  12,  15,  22,  and  26  are  the  centers  for  these  curves. 

Each  point  input  defines  the  point  number,  followed  by  its  x  and  y  position 
for  a  two-dimension  plane  strain  model,  or  its  z  and  r  position,  if  the 
model  is  axisymmetric.  For  this  example,  the  projectile  is  an 
axisymmetric  structure,  so  the  coordinates  are  in  terms  of  the 
longitudinal  position  and  the  radial  position  of  the  point. 
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Figure  3.0 

Geometry  Sub-Area  Definitions 


Following  the  corner  point  definitions,  the  perimeter  lines  for  each 
of  the  sub-areas  are  defined.  Table  3.2  gives  the  line  inputs  for  this 
example,  and  Figure  3.2  shows  the  line  definitions  graphically.  Each  line 
input  defines  the  line  number,  followed  by  the  two  line  end  points,  the 
number  of  finite  divisions  for  this  line,  and  the  center  point  if  the  line  is 
a  radius.  When  defining  the  lines,  as  when  defining  the  sub-areas  and 
their  corner  points,  the  modeler  needs  to  have  a  scheme  in  mind  for  how 
the  mesh  is  to  be  generated.  Line  divisions  should  be  consistent  for 
opposite  sides  of  each  sub-area,  so  that  a  consistent  quadrilateral  grid 
will  be  generated.  In  addition,  the  line  divisions  for  adjacent  areas  should 
be  consistent  so  that  the  entire  structure  can  be  merged  together. 
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Table  3.1 

Sub-Area  Corner  Point  Inputs 

P.1.0.0.0 
P.2.0.0.5 
P, 3. 2.0,0. 0 
P,4.2.0,0.5 
P.5.4.25,0.0 
P.6,4.25,0.5 
P, 7,4. 25,1.0 
P.8,2.0,0.65 
P.9,5-7i.0.0 
P. 10,5. 7i, 0.5 
P.11,4.5,1.25 
P.12,4.0443,1.4557 
P.13,4.5,1.5 
P.14,5.5,1.5 
P. 15,5.9557,1.4557 
P.16.5.5,1.25 
P, 17,5.75,1.0 
P, 18,8.0,6.0 
P, 19,8.0,0.5 
P.20.8.0,0.7 
P.21.8.5,0.825 
P, 22,8.0600,1. 7 
P.23.8.5. 0.6 
P.24,8.5,0.5 
P, 25,8.75,0.75 
P, 26.8.9557,0.2943 
P, 27,9.6,1.35 
P, 28,9.5  1.5 
P, 29,10.6,1.35 
P-,30,10.0,1.5 
P, 31. 10.5,0.0 
P.32,10.5,0.5 

Figure  3. 1 

Points  Defining  the  Sub-Areas  of  the  Structure 
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One  may  have  noticed  that  three  of  the  sub-areas  are  defined  by 
more  than  four  lines  and  more  than  four  perimeter  points.  Since  typical 
structures  involve  complicated  shapes,  which  are  not  easily  described  by 
four  nearly  parallel  sides,  mesh  generation  requires  the  use  of  splines. 
Splines  are  groups  of  lines,  which  are  to  be  treated  as  a  single  line  in  the 
mesh  generation  process.  This  approach  makes  mesh  generation  of 
complicated  shapes  easier  to  accomplish.  Table  3.3  shows  the  five  spline 
inputs  for  this  example.  When  defining  splines,  at  most  five  lines  can  be 
connected,  and  they  should  be  listed  in  the  order  they  appear  in  the 
geometry.  As  with  the  definition  of  individual  lines,  the  total  number  of 
divisions  for  a  spline  should  be  consistent  with  the  opposite  side  of  the 
sub-area  which  it  describes. 


Table  3.3 
Spline  Inputs 


5.1.21.23.24.. , 

5.2.28.27.26.. , 

5.3.32.33.. ., 

5.4.34.35.. ., 

5.5.37.38.. ., 


Following  the  spline  inputs,  the  model  is  ready  to  be  meshed  with  a 
finite  element  grid.  The  mesh  commands  define  each  sub-area  to  be 
meshed  in  accordance  with  the  points,  lines,  and  splines  described  earlier. 
For  each  area,  the  four  corner  points  are  entered  in  either  a  clock-wise  or 
counter  clock-wise  direction,  depending  on  the  convention  used  in  the 
finite  element  software.  Table  3.4  shows  the  area  mesh  commands  used 
in  this  example.  The  command  line  consists  of  the  area  number,  followed 
by  the  four  corner  points,  followed  by  the  mesh  technique  (1  or  2),  and 
finally  the  material  number,  so  that  material  properties  can  be 
discriminated  later  in  the  finite  element  analysis.  Two  mesh  technique 
numbers  were  developed  in  the  software,  since  any  one  technique  is  not 
always  the  best,  depending  on  the  shape  of  the  sub-area.  Technique  1 
works  well  for  areas  which  are  more  square  or  rectangular,  as  opposed  to 
areas  which  include  a  curved  perimeter.  Some  experimentation  is  required 
to  develop  the  best  looking  grid  for  a  particular  sub-area.  The  algorithm 
is  adequately  general  to  permit  the  grid  to  take  several  forms.  Results 
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will  vary  depending  on  the  direction  and  the  order  used  for  the  area  corner 
points  in  the  mesh  command.  The  modeler  needs  to  decide  which  results 
are  best  for  the  particular  application.  Figure  3.3  shows  the  resulting 
mesh  using  the  commands  in  Table  3.4.  Figure  3.4  shows  a  close-up  of  the 
mesh  in  the  bulkhead  region.  Figure  3.5  shows  a  close-up  of  the  root  of 
the  front  scoop  with  a  point  and  line  plot.  Figure  3.6  is  the  resulting  mesh 
for  the  front  scoop.  The  mesh  generation  for  both  of  these  areas  included 
the  use  of  splines. 


Table  3.4 

Area  Mesh  Command  Inputs 


Figure  3.3 

The  Complete  Finite  Element  Grid 


Figure  3.4 

Close-Up  of  the  Bulkhead  Region  Grid 
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Figure  3.5 

Point  and  Line  Piot  of  the  Front  Scoop  Region 

22 
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Figure  3.6 

Front  Scoop  Region  Grid 
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4.  DETAILS  OF  THE  MESH  GENERATION  PROCESS 


For  each  sub-area  mesh  command,  the  mesh  generation  process 
begins  by  identifying  the  line  or  spline  which  connects  the  first  two 
points  in  the  command  line.  Since  the  line  or  spline  could  have  been 
defined  in  the  opposite  direction,  a  test  is  performed  on  each  line  to  see  if 
it  is  reversed.  If  the  two  points  define  a  spline,  each  line  in  the  spline  is 
identified  as  well.  The  meshing  process  begins  by  subdividing  all  the  edge 
lines  or  splines  according  to  the  number  of  divisions  in  the  line 
definitions.  Several  scratch  files  are  also  maintained  to  keep  track  of  the 
line  division  node  numbers  for  future  reference,  particularly  when 
merging  adjacent  areas.  The  nodes  defined  on  the  first  line  of  the  area  are 
then  swept  across  the  width  of  the  area  based  on  the  division  spacing  of 
the  line  or  spline  between  the  second  and  third  points.  This  where  the 
algorithm  uses  either  technique  1  or  2  when  generating  the  grid.  If 
technique  1  is  used,  the  nodes  sweep  straight  across  to  the  opposite  side. 
With  technique  2,  the  curvature  of  the  second  side  adjusts  the  node 
spacing  as  they  sweep  across.  Within  each  sub-area  all  elements  are 
connected  by  adjacent  nodes  according  to  the  number  of  divisions  on  the 
perpendicular  sides,  but  adjacent  areas  are  not  yet  connected.  Following 
completion  of  the  meshing  of  all  sub-areas,  the  operator  can  merge 
perimeter  nodes,  so  that  a  solid  structure  is  created.  Merging  is 
accomplished  by  cycling  through  all  sub-areas  to  identify  those  which 
happened  to  have  common  boundary  lines.  Scratch  files,  created 
previously,  contain  the  edge  nodes  for  each  line.  The  merge  routine 
selects  the  lowest  node  number  and  makes  it  common  to  all  elements 
connecting  at  that  location.  After  merging,  some  node  numbers  will  no 
longer  be  in  use.  The  operator  can  then  select  the  compress  command, 
which  will  identify  ail  un-used  nodes  in  the  node  array.  Un-used  nodes  are 
then  re-assigned  to  elements  having  higher  node  definitions.  This 
command  is  useful  in  ensuring  that  the  size  of  the  stiffness  matrix  is 
minimized  in  the  finite  element  analysis.  Figure  4.1  shows  a  close-up  of 
the  point  and  line  plot  of  the  left  end  of  the  example  structure.  Figure  4.2 
shows  the  node  and  element  definitions  at  the  vertical  sub-area  boundary. 
At  this  boundary,  the  nodes  have  been  merged  so  that  the  elements  are 
structurally  connected. 
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Figure  4. 1 

Close-Up  of  the  Left  End  of  the  Example  Structure 


Figure  4.2 

Merged  Node  and  Element  Plot  at  the  Sub-Area  Boundary 
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APPENDIX  A 

EXAMPLE  INPUT  AND  OUTPUT 


A.1  -  INPUT  FILE 

A.2  -  NODE  FILE  OUTPUT 

A.3  ••  ELEMENT  FILE  OUTPUT 
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A.1  -  INPUT  FILE 


P.2.0,0.5 
P, 3, 2. 0,0.0 
P, 4, 2. 0,0. 5 
P,5,4.25,0.0 
P,6,4.25,0.5 
P, 7,4. 25.1.0 
P,8,2.q.0.65 
P, 9.5. 7^0.0 
P, 10,5.75,0.5 
P, 11, 4. 5, 1.25 
P, 12,4. 0443, 1.4557 
P, 13,4. 5, 1.5 
P, 14,5. 5,1. 5 
P, 15, 5.9557,1. 4557 
P,16,5.^1.25 
P,17,5. 75,1.0 
P, 18,8.0,0.0 
P,19,8.0,0.5 
P,20,8.0,0.7 
P,21. 8.5,0.825 
P,22,8.0000,1.7 
P, 23,8.5, 0.6 
P,24.8.5,0.5 
P, 25, 8.75,0. 75 
P, 26,8.9557,0. 2943 
P, 27,9.6, 1.35 
P.28,9.5,1.5 
P,29,10.6.1.35 
P.30.10.0,1.5 
P.31.10.5,0.0 
P, 32,10.5, 0.5 

•-.2. 3, 4, 5,, 

4. 3. 4. 2. 5. , 

L. 4. 1.2,5., 

L, 5. 3.5, 20,. 

1.. 6.5.6.5., 

1.7.6.4.26., 

L,8.5.9  20.. 

L, 9.9,16,5,. 

1.16.10.6.66., 

L,11.9,!8,20 

l. 12.16,16,5, 

1.13.19.10.26., 

L. 14, 18,23,8,, 

L. 15,23.24,5.. 

1.16.24.19.8., 

L. 17,23,31,4,, 

1.18.31.32.5., 

1.20.4.8.5., 

L, 21.6,7.5  , 

1.22.7.8.26., 

L,23,7,l1.4j2 

L,24,l1,l5.S.. 

4.25.13.14.26., 

1.26.14.16.3., 
1,27,16,17,4,15 
L.28,10,17,5,, 

4.29.19.20.5., 

4.30.20.17.26., 

4,31,24,25,5,66 

4.32.25.21.3., 
4,33,21,20,5.22 

4.34.25.27.16., 

4.35.27.29.3., 

4.36.29.30.3., 

4.37.30.28.3., 

4.38.28.21.16., 

5.1.61.63.64.. , 

5.2.28.27.26.. , 

5.3.32.33.. ., 
S/^f34,35, , , , 

5.5.37.38.. ., 
A,1,3.4,2,1,<,1 
A,2,3,4,6,5,1,1 
A,3,5,9. 16.6, 1,1 
a.4,9.i6.i6,1oj,i 
A,5.l6.2i,24.l6,1,1 
A,6,23, 31,32, 24, 1,1 
A, 7,4,6. 7, 8. 1.2 
A.8,6,i6J4.i6.1,2 
A,9,i6,i6.26,17,1,2 
A.16.26.2S, 24. 16,6,2 
A.11, 29,30,21, 25, 1.2 


A.2  -  NODE  FILE  OUTPUT 


1 

2.00000 

0.000000 

94 

3.01250 

0.300000 

2 

2.00000 

0.100000 

95 

3.01250 

0.400000 

3 

2.00000 

0.200000 

96 

3.01250 

0.500000 

4 

2.00000 

0.300000 

97 

3.12500 

0.000000 

S 

2.00000 

0.400000 

98 

3.12500 

0.100000 

6 

2.00000 

0.500000 

99 

3.12500 

0.200000 

7 

1.60000 

0.000000 

100 

3.12500 

0.300000 

8 

1.60000 

0.100000 

101 

3.12500 

0.400000 

9 

1.60000 

0.200000 

102 

3.12500 

0.500000 

10 

1.60000 

0.300000 

103 

3.23750 

0.000000 

11 

1.60000 

0.400000 

104 

3.23750 

0.100000 

12 

1.60000 

0.500000 

105 

3.23750 

0.200000 

13 

1.20000 

0.000000 

106 

3.23750 

0.300000 

14 

1.20000 

0.100000 

107 

3.23750 

0.400000 

15 

1.20000 

0.200000 

108 

3.23750 

0.500000 

16 

1.20000 

0.300000 

109 

3.35000 

0.000000 

17 

1.20000 

0.400000 

110 

3.35000 

0.100000 

18 

1.20000 

0.500000 

111 

3.35000 

0.200000 

19 

0.800000 

0.000000 

112 

3.35000 

0.300000 

20 

0.800000 

0.100000 

113 

3.35000 

0.400000 

21 

0.800000 

0.200000 

114 

3.35000 

0.500000 

22 

0.800000 

0.300000 

115 

3.46250 

0.000000 

23 

0.800000 

0.400000 

116 

3.46250 

0.100000 

24 

0,800000 

0.500000 

117 

3.46250 

0.200000 

25 

0.400000 

0.000000 

118 

3.46250 

0.300000 

26 

0.400000 

0.100000 

119 

3.46250 

0.400000 

27 

0.400000 

0.200000 

120 

3.46250 

0.500000 

28 

0.400000 

0.300000 

121 

3.57500 

0.000000 

29 

0.400000 

0.400000 

122 

3.57500 

0.100000 

30 

0.400000 

0.500000 

123 

3.57500 

0.200000 

31 

0.000000 

0.000000 

124 

3.57500 

0.300000 

32 

0.000000 

0.100000 

125 

3.57500 

0.400000 

33 

0.000000 

0.200000 

126 

3.57500 

0.500000 

34 

0,000000 

0.300000 

127 

3.68750 

0.000000 

35 

0.000000 

0.400000 

128 

3.68750 

0.100000 

36 

0.000000 

o.sooooo 

129 

3.68750 

0.200000 

37 

2.00000 

0.000000 

130 

3.68750 

0.300000 

38 

2.00000 

0.100000 

131 

3.68750 

0.400000 

39 

2.00000 

0.200000 

132 

3.68750 

0.500000 

40 

2.00000 

0.300000 

133 

3.80000 

0.000000 

41 

2.00000 

0.400000 

134 

3.80000 

0.100000 

42 

2.00000 

0.500000 

135 

3.80000 

0.200000 

43 

2.11250 

0.000000 

136 

3.80000 

0.300000 

44 

2.11250 

0.100000 

137 

3.80000 

0.400000 

45 

2.11250 

0.200000 

138 

3.80000 

0.500000 

46 

2.11250 

0.300000 

139 

3.91250 

0.000000 

47 

2.11250 

0.400000 

140 

3.91250 

0.100000 

48 

2,11250 

0.500000 

141 

3.91250 

0.200000 

49 

2.22500 

0.000000 

142 

3.91250 

0.300000 

50 

2.22500 

0.100000 

143 

3.91250 

0.400000 

51 

2.22500 

0.200000 

144 

3.91250 

0.500000 

52 

2.22500 

0.300000 

145 

4.02500 

0.000000 

53 

2.22500 

0.400000 

146 

4.02500 

0.100000 

54 

2.22500 

0.500000 

147 

4.02500 

0.200000 

55 

2.33750 

0.000000 

148 

4.02500 

0.300000 

56 

2.33750 

0.100000 

149 

4.02500 

0.400000 

57 

2.33750 

0.200000 

150 

4.02500 

0.500000 

58 

2.33750 

0.300000 

151 

4.13750 

0.000000 

59 

2.33750 

0.400000 

152 

4.13750 

0.100000 

60 

2.33750 

0.500000 

153 

4.13750 

0.200000 

61 

2.45000 

0.000000 

154 

4.13750 

0.300000 

62 

2.45000 

0.100000 

155 

4.13750 

0.400000 

63 

2.45000 

0.200000 

156 

4.13750 

0.500000 

64 

2.45000 

0.300000 

157 

4.25000 

0.000000 

65 

2.45000 

0.400000 

158 

4.25000 

0.100000 

66 

2.45000 

0.500000 

159 

4.25000 

0.200000 

67 

2.56250 

0.000000 

160 

4.25000 

0.300000 

68 

2.56250 

0.100000 

161 

4.25000 

0.400000 

69 

2.56250 

0.200000 

162 

4.25000 

0.500000 

70 

2.56250 

0.300000 

163 

4.25000 

0.000000 

71 

2.56250 

0.400000 

164 

4.32500 

0.000000 

72 

2.56250 

0.500000 

165 

4.40000 

0.000000 

73 

2.67500 

0.000000 

166 

4.47500 

0.000000 

74 

2.67500 

0.100000 

167 

4.55000 

0.000000 

75 

2.67500 

0.200000 

168 

4.62500 

0.000000 

76 

2.67500 

0.300000 

169 

4.70000 

0.000000 

77 

2.67500 

0.400000 

170 

4.77500 

0.000000 

78 

2.67500 

0.500000 

171 

4.85000 

0.000000 

79 

2.78750 

0.000000 

172 

4.92500 

0.000000 

80 

2.78750 

0.100000 

173 

5.00000 

0.000000 

81 

2.78750 

0.200000 

174 

5.07500 

O.MOOOO 

82 

2.78750 

0.300000 

175 

5.15000 

0.000000 

83 

2.78750 

0.400000 

176 

5.22500 

0.000000 

84 

2.78750 

0.500000 

177 

5.30000 

0.000000 

85 

2.90000 

0.000000 

178 

5.37500 

0.000000 

86 

2.90000 

0.100000 

179 

5.45000 

0.000000 

87 

2.90000 

0.200000 

180 

5.52500 

0.000000 

88 

2.90000 

0.300000 

181 

5.60000 

0.000000 

89 

2.90000 

0.400000 

182 

5.67500 

0.000000 

90 

2.90000 

0.500000 

183 

5.75000 

0.000000 

91 

3.01250 

0.000000 

184 

4.25000 

0.100000 

92 

3.01250 

0.100000 

185 

4.32500 

0.100000 

93 

3.01250 

0.200000 

186 

4.40000 

0.100000 

187 

4.47500 

0.100000 

280 

5.15000 

0.500000 

188 

4.55000 

0.100000 

281 

5.22500 

0.500000 

189 

4.62500 

0.100000 

Z6Z 

5.30000 

0.500000 

190 

4.70000 

0.100000 

283 

5.37500 

0.500000 

191 

4.77500 

0.100000 

284 

5.45000 

0.500000 

192 

4.85000 

0.100000 

285 

5.52500 

0.500000 

193 

4.92500 

0.100000 

286 

5.60000 

0.500000 

194 

5.00000 

0.100000 

287 

5.67500 

0.500000 

195 

5.07500 

0.100000 

288 

5.75000 

0.500000 

196 

5.15000 

0.100000 

289 

5.75000 

0.000000 

197 

5.22500 

0.100000 

290 

5.86250 

0.000000 

198 

5.30000 

0.100000 

291 

5.97500 

0.000000 

199 

5.37500 

0.100000 

292 

6.08750 

0.000000 

200 

5.45000 

0.100000 

293 

6.20000 

0.000000 

201 

5.52500 

0.100000 

294 

6.31250 

0.000000 

202 

5.60000 

0.100000 

295 

6.42500 

0.000000 

203 

5.67500 

0.100000 

296 

6.53750 

0.000000 

204 

5.75000 

0.100000 

297 

6.65000 

0.000000 

205 

4.25000 

0.200000 

298 

6.76250 

0.000000 

206 

4.32500 

0.200000 

299 

6.87500 

0.000000 

207 

4.40000 

0.200000 

300 

6.98750 

0.000000 

208 

4.47500 

0.200000 

301 

7.10000 

0.000000 

209 

4.55000 

0.200000 

302 

7.21250 

0.000000 

210 

4.62500 

0.200000 

303 

7.32500 

0.000000 

211 

4.70000 

0.200000 

304 

7.43750 

0.000000 

212 

4.77500 

0.200000 

305 

7.55000 

0.000000 

213 

4.85000 

0.200000 

306 

7.66250 

0.000000 

214 

4.92500 

0.200000 

307 

7.77500 

0.000000 

215 

5.00000 

0.200000 

308 

7.88750 

0.000000 

216 

5.07500 

0.200000 

309 

8.00000 

0.000000 

217 

5.15000 

0.200000 

310 

5.75000 

0.100000 

218 

5.22500 

0.200000 

311 

5.86250 

0.100000 

219 

5.30000 

0.200000 

312 

5.97500 

0.100000 

220 

5.37500 

0.200000 

313 

6.08750 

0.100000 

221 

5.45000 

0.200000 

314 

6.20000 

0.100000 

222 

5.52500 

0.200000 

315 

6.31250 

0.100000 

223 

5.60000 

0.200000 

316 

6.42500 

0.100000 

224 

5.67500 

0.200000 

317 

6.53750 

0.100000 

225 

5.75000 

0.200000 

318 

6.65000 

0.100000 

226 

4.25000 

0.300000 

319 

6.76250 

0.100000 

227 

4.32500 

0.300000 

320 

6.87500 

0.100000 

228 

4.40000 

0.300000 

321 

6.98750 

0.100000 

229 

4.47500 

0.300000 

322 

7.10000 

0.100000 

230 

4.55000 

0.300000 

323 

7.21250 

0.100000 

231 

4.62500 

0.300000 

324 

7.32500 

0.100000 

232 

4.70000 

0.300000 

325 

7.43750 

0.100000 

233 

4.77500 

0.300000 

326 

7.55000 

0.100000 

234 

4.85000 

0.300000 

327 

7.66250 

0.100000 

235 

4.92500 

0.300000 

328 

7.77500 

0.100000 

236 

5.00000 

0.300000 

329 

7.88750 

0.100000 

237 

5.07500 

0.300000 

330 

8.00000 

0.100000 

238 

5.15000 

0.300000 

331 

5.75000 

0.200000 

239 

5.22500 

0.300000 

332 

5.86250 

0.200000 

240 

5.30000 

0.300000 

333 

5.97500 

0.200000 

241 

5.37500 

0.300000 

334 

6.08750 

0.200000 

242 

5.45000 

0.300000 

335 

6.20000 

0.200000 

243 

5.52500 

0.300000 

336 

6.31250 

0.200000 

244 

5.60000 

0.300000 

337 

6.42500 

0.200000 

245 

5.67500 

0.300000 

338 

6.53750 

0.200000 

246 

5.75000 

0.300000 

339 

6.65000 

0.200000 

247 

4.25000 

0.400000 

340 

6.76250 

0.200000 

248 

4.32500 

0.400000 

341 

6.87500 

0.200000 

249 

4.40000 

0.400000 

342 

6.98750 

0.200000 

250 

4.47500 

0.400000 

343 

7.10000 

0.200000 

251 

4.55000 

0.400000 

344 

7.21250 

0.200000 

252 

4.62500 

0.400000 

345 

7.32500 

0.200000 

253 

4.70000 

0.400000 

346 

7.43750 

0.200000 

254 

4.77500 

0.400000 

347 

7.55000 

0.200000 

255 

4.85000 

0.400000 

348 

7.66250 

0.200000 

256 

4.92500 

0.400000 

349 

7.77500 

0.200000 

257 

5.00000 

0.400000 

350 

7.88750 

0.200000 

258 

5.07500 

0.400000 

351 

8.00000 

0.200000 

259 

5.15000 

0.400000 

352 

5.75000 

0.300000 

260 

5.22500 

0.400000 

353 

5.86250 

0.300000 

261 

5.30000 

0.400000 

354 

5.97500 

0.300000 

262 

5.37500 

0.400000 

355 

6.08750 

0.300000 

263 

5.45000 

0.400000 

356 

6.20000 

0.300000 

264 

5.52500 

0.400000 

357 

6.31250 

0.300000 

265 

5.60000 

0.400000 

358 

6.42500 

0.300000 

266 

5.67500 

0.400000 

359 

6.53750 

0.300000 

267 

5.75000 

0.400000 

360 

6.65000 

0.300000 

268 

4.25000 

0.500000 

361 

6.76250 

0.300000 

269 

4.32500 

0.500000 

362 

6.87500 

0.300000 

270 

4.40000 

0.500000 

363 

6.98750 

0.300000 

271 

4.47500 

0.500000 

364 

7.10000 

0.300000 

272 

4.55000 

0.500000 

365 

7.21250 

0.300000 

273 

4.62500 

0.500000 

366 

7.32500 

0.300000 

274 

4.70000 

0.500000 

367 

7.43750 

0.300000 

275 

4.77500 

0.500000 

368 

7.55000 

0.300000 

276 

4.85000 

0.500000 

369 

7.66250 

0.300000 

277 

4.92500 

0.500000 

370 

7.77500 

0.300000 

278 

5.00000 

0.500000 

371 

7.88750 

0.300000 

279 

5.07500 

0.500000 

372 

8.00000 

0.300000 

373 

5.75000 

0.400000 

466 

8.37500 

0.500000 

37i 

5.86250 

0.400000 

467 

8.43750 

0.500000 

375 

5.97500 

0.400000 

468 

8.50000 

0.500000 

376 

6.08750 

0.400000 

469 

8.50000 

0.000000 

377 

6.20000 

0.400000 

470 

9.00000 

0.000000 

378 

6.31250 

0.400000 

471 

9.50000 

0.000000 

379 

6.42500 

0.400000 

472 

10.0000 

0.000000 

380 

6.53750 

0.400000 

473 

10.5000 

0.000000 

381 

6.65000 

0.400000 

474 

8.50000 

0.100000 

382 

6.76250 

0.400000 

475 

9.00000 

0.100000 

383 

6.87500 

0.400000 

476 

9.50000 

0.100000 

384 

6.98750 

0.400000 

477 

10.0000 

0.100000 

385 

7.10000 

0.400000 

478 

10.5000 

0.100000 

386 

7.21250 

0.400000 

479 

8.50000 

0.200000 

387 

7.32500 

0.400000 

480 

9.00000 

0.200000 

388 

7.43750 

0.400000 

481 

9.50000 

0.200000 

389 

7.55000 

0.400000 

482 

10.0000 

0.200000 

390 

7.66250 

0.400000 

483 

10.5000 

0.200000 

391 

7.77500 

0.400000 

484 

8.50000 

0.300000 

392 

7.88750 

0.400000 

485 

9.00000 

0.300000 

393 

8.00000 

0.400000 

486 

9.50000 

0.300000 

394 

5.75000 

0.500000 

487 

10.0000 

0.300000 

395 

5.86250 

0.500000 

488 

10.5000 

0.300000 

396 

5.97500 

0.500000 

489 

8.50000 

0.400000 

397 

6.08750 

0.500000 

490 

9.00000 

0.400000 

398 

6.20000 

0.500000 

491 

9.50000 

0.400000 

399 

6.31250 

0.500000 

492 

10.0000 

0.400000 

400 

6.42500 

0.500000 

493 

10.5000 

0.400000 

401 

6.53750 

0.500000 

494 

8.50000 

0.500000 

402 

6.65000 

0.500000 

495 

9.00000 

0.500000 

403 

6.76250 

0.500000 

496 

9.50000 

0.500000 

404 

6.87500 

0.500000 

497 

10.0000 

0.500000 

405 

6.98750 

0.500000 

498 

10.5000 

0.500000 

406 

7.10000 

0.500000 

499 

2.00000 

0.500000 

407 

7.21250 

0.500000 

500 

2.11250 

0.500000 

408 

7.32500 

0.500000 

501 

2.22500 

0.500000 

409 

7.43750 

0.500000 

502 

2.33750 

0.500000 

410 

7.55000 

0.500000 

503 

2.45000 

0.500000 

411 

7.66250 

0.500000 

504 

2.56250 

0.500000 

412 

7.77500 

0.500000 

505 

2.67500 

0.500000 

413 

7.88750 

0.500000 

506 

2.78750 

0.500000 

414 

8.00000 

0.500000 

507 

2.90000 

0.500000 

415 

8.00000 

0.000000 

508 

3.01250 

0.500000 

416 

8.06250 

0.000000 

509 

3.12500 

0.500000 

417 

8.12500 

0.000000 

510 

3.23750 

0.500000 

418 

8.18750 

0.000000 

511 

3.35000 

0.500000 

419 

8.25000 

0.000000 

512 

3.46250 

0.500000 

420 

8.31250 

0.000000 

513 

3.57500 

0.500000 

421 

8.37500 

0.000000 

514 

3.68750 

0.500000 

422 

8.43750 

0.000000 

515 

3.80000 

0.500000 

423 

8.50000 

0.000000 

516 

3.91250 

0.500000 

424 

8.00000 

0.100000 

517 

4.02500 

0.500000 

425 

8.06250 

0.100000 

518 

4.13750 

0.500000 

426 

8.12500 

0.100000 

519 

4.25000 

0.500000 

427 

8.18750 

0.100000 

520 

2.00000 

0.530000 

428 

8.25000 

0.100000 

521 

2.11250 

0.533500 

429 

8.31250 

0.100000 

522 

2.22500 

0.537000 

430 

8.37500 

0.100000 

523 

2.33750 

0.540500 

431 

8.43750 

0.100000 

524 

2.45000 

0.544000 

432 

8.50000 

0.100000 

525 

2.56250 

0.547500 

433 

8.00000 

0.200000 

526 

2.67500 

0.551000 

434 

8.06250 

0.200000 

527 

2.78750 

0.554500 

435 

8.12500 

0.200000 

528 

2.90000 

0.558000 

436 

8.18750 

0.200000 

529 

3.01250 

0.561500 

437 

8.25000 

0.200000 

530 

3.12500 

0.565000 

438 

8.31250 

0.200000 

531 

3.23750 

0.568500 

439 

8.37500 

0.200000 

532 

3.35000 

0.572000 

440 

8.43750 

0.200000 

533 

3.46250 

0.575500 

441 

8.50000 

0.200000 

534 

3.57500 

0.579000 

442 

8.00000 

0.300000 

535 

3.68750 

0.582500 

443 

8.06250 

0.300000 

536 

3.80000 

0.586000 

444 

8.12500 

0.300000 

537 

3.91250 

0.589500 

445 

8.18750 

0.300000 

538 

4.02500 

0.593000 

446 

8.25000 

0.300000 

539 

4.13750 

0.596500 

447 

8.31250 

0.300000 

540 

4.25000 

0.600000 

448 

8.37500 

0.300000 

541 

2.00000 

0.560000 

449 

8.43750 

0.300000 

542 

2.11250 

0.567000 

450 

8.50000 

0.300000 

543 

2.22500 

0.574000 

451 

8.00000 

0.400000 

544 

2.33750 

0.581000 

452 

8.06250 

0.400000 

545 

2.45000 

0.588000 

453 

8.12500 

0.400000 

546 

2.56250 

0.595000 

454 

8.18750 

0.400000 

547 

2.67500 

0.602000 

455 

8.25000 

0.400000 

548 

2.78750 

0.609000 

456 

8.31250 

0.400000 

549 

2.90000 

0.616000 

457 

8.37500 

0.400000 

550 

3.01250 

0.623000 

458 

8.43750 

0.400000 

551 

3.12500 

0.630000 

459 

8.50000 

0.400000 

552 

3.23750 

0.637000 

460 

8.00000 

0.500000 

553 

3.35000 

0.644000 

461 

8.06250 

0.500000 

554 

3.46250 

0.651000 

462 

8.12500 

0.500000 

555 

3.57500 

0.658000 

463 

8.18750 

0.500000 

556 

3.68750 

0.665000 

464 

8.25000 

0.500000 

557 

3.80000 

0.672000 

465 

8.31250 

0.500000 

558 

3.91250 

0.679000 

559 

4.02500 

0.686000 

652 

4.70000 

0.600000 

560 

4.13750 

0.693000 

653 

4.77500 

0.600000 

561 

4.25000 

0.700000 

654 

4.85000 

0.600000 

562 

2.00000 

0.590000 

655 

4.92500 

0.600000 

563 

2.11250 

0.600500 

656 

5.00000 

0.600000 

564 

2.22500 

0.611000 

657 

5.07500 

0.600000 

565 

2.33750 

0.621500 

658 

5.15000 

0.600000 

566 

2.45000 

0.632000 

659 

5.22500 

0.600000 

567 

2.56250 

0.642500 

660 

5.30000 

0.600000 

568 

2.67500 

0.653000 

661 

5.37500 

0.600000 

569 

2.78750 

0.663500 

662 

5.45000 

0.600000 

570 

2.90000 

0.674000 

663 

5.52500 

0.600000 

571 

3.01250 

0.684500 

664 

5.60000 

0.600000 

572 

3.12500 

0.695000 

665 

5.67500 

0.600000 

573 

3.23750 

0.705500 

666 

5.75000 

0.600000 

574 

3.35000 

0.716000 

667 

4.25000 

0.700000 

575 

3.46250 

0.726500 

668 

4.32500 

0.700000 

576 

3.57500 

0.737000 

669 

4.40000 

0.700000 

577 

3.68750 

0.747500 

670 

4.47500 

0.700000 

578 

3.80000 

0.758000 

671 

4.55000 

0.700000 

579 

3.91250 

0.768500 

672 

4.62500 

0.700000 

580 

4.02500 

0.779000 

673 

4.70000 

0.700000 

581 

4.13750 

0.789500 

674 

4.77500 

0.700000 

582 

4.25000 

0.800000 

675 

4.85000 

0.700000 

583 

2.00000 

0.620000 

676 

4.92500 

0.700000 

584 

2.11250 

0.634000 

677 

5.00000 

0.700000 

585 

2.22500 

0.648000 

678 

5.07500 

0.700000 

586 

2.33750 

0.662000 

679 

5.15000 

0.700000 

587 

2.45000 

0.676000 

680 

5.22500 

0.700000 

588 

2.56250 

0.690000 

681 

5.30000 

0.700000 

589 

2.67500 

0.704000 

682 

5.37500 

0.700000 

590 

2.78750 

0.718000 

683 

5.45000 

0.700000 

591 

2.90000 

0.732000 

684 

5.52500 

0.700000 

592 

3.01250 

0.746000 

685 

5.60000 

0.700000 

593 

3.12500 

0.760000 

686 

5.67500 

0.700000 

594 

3.23750 

0.774000 

687 

5.75000 

0.700000 

595 

3.35000 

0.788000 

688 

4.25000 

0.800000 

596 

3.46250 

0.802000 

689 

4.32500 

0.800000 

597 

3.57500 

0.816000 

690 

4.40000 

0.800000 

598 

3.68750 

0.830000 

691 

4.47500 

0.800000 

599 

3.80000 

0.844000 

692 

4.55000 

0.800000 

600 

3.91250 

0.858000 

693 

4.62500 

0.800000 

601 

4.02500 

0.872000 

694 

4.70000 

0.800000 

602 

4.13750 

0.886000 

695 

4.77500 

0.800000 

603 

4.25000 

0.900000 

696 

4.85000 

0.800000 

604 

2.00000 

0.650000 

697 

4.92500 

0.800000 

60S 

2.11250 

0.667500 

698 

5.00000 

0.800000 

606 

2.22500 

0.685000 

699 

5.07500 

0.800000 

607 

2.33750 

0.702500 

700 

5.15000 

0.800000 

608 

2.45000 

0.720000 

701 

5.22500 

0.800000 

609 

2.56250 

0.737500 

702 

5.30000 

0.800000 

610 

2.67500 

0.755000 

703 

5.37500 

0.800000 

611 

2.78750 

0.772500 

704 

5.45000 

0.800000 

612 

2.90000 

0.790000 

705 

5.52500 

0.800000 

613 

3.01250 

0.807500 

706 

5.60000 

0.800000 

614 

3.12500 

0.825000 

707 

5.67500 

0.800000 

615 

3.23750 

0.842500 

708 

5.75000 

0.800000 

616 

3.35000 

0.860000 

709 

4.25000 

0.900000 

617 

3.46250 

0.877500 

710 

4.32500 

0.900000 

618 

3.57500 

0.895000 

711 

4.40000 

0.900000 

619 

3.68750 

0.912500 

712 

4.47500 

0.900000 

620 

3.80000 

0.930000 

713 

4.55000 

0.900000 

621 

3.91250 

0.947500 

714 

4.62500 

0.900000 

622 

4.02500 

0.965000 

715 

4.70000 

0.900000 

623 

4.13750 

0.982500 

716 

4.77500 

0.900000 

624 

4.25000 

1 .00000 

717 

4.85000 

0.900000 

625 

4.25000 

0.500000 

718 

4.92500 

0.900000 

626 

4.32500 

0.500000 

719 

5.00000 

0.9(K)000 

627 

4.40000 

0.500000 

720 

5.07500 

0.900000 

628 

4.47500 

0.500000 

721 

5.15000 

0.900000 

629 

4.55000 

0.500000 

722 

5.22500 

0.900000 

630 

4.62500 

0.500000 

723 

5.30000 

0.900000 

631 

4.70000 

0.500000 

724 

5.37500 

0.900000 

632 

4.77500 

0.500000 

725 

5.45000 

0.900000 

633 

4.85000 

0.500000 

726 

5.52500 

0.900000 

634 

4.92500 

0.500000 

727 

5.60000 

0.900000 

635 

5.00000 

0.500000 

728 

5.67500 

0.900000 

636 

5.07500 

0.500000 

729 

5.75000 

0.900000 

637 

5.15000 

0.500000 

730 

4.25000 

0.999999 

638 

5.22500 

0.500000 

731 

4.32500 

0.999999 

639 

5.30000 

0.500000 

732 

4.40000 

0.999999 

640 

5.37500 

0.500000 

733 

4.47500 

0.999999 

641 

5.45000 

0.500000 

734 

4.55000 

0.999999 

642 

5.52500 

0.500000 

735 

4.62500 

0.999999 

643 

5.60000 

0.500000 

736 

4.70000 

1.00000 

644 

5.67500 

0.500000 

737 

4.77500 

1.00000 

645 

5.75000 

0.500000 

738 

4.85000 

1.00000 

646 

4.25000 

0.600000 

739 

4.92500 

1.00000 

647 

4.32500 

0.600000 

740 

5.00000 

I.OOIKK) 

648 

4.40000 

0.600000 

741 

5.07500 

1.00000 

649 

4.47500 

0.600000 

742 

5.15000 

1.00000 

650 

4.55000 

0.60(»00 

743 

5.22500 

1.00000 

651 

4.62500 

0.600000 

744 

5.30000 

1.00000 

745 

5.37500 

1.00000 

838 

4.65000 

1.33333 

746 

5.45000 

1.00000 

839 

4.70000 

1 .33333 

747 

5.52500 

1.00000 

840 

4.75000 

1 .33333 

74« 

5.60000 

1.00000 

841 

4.80000 

1.33333 

749 

5.67500 

1.00000 

842 

4.85000 

1 .33333 

750 

5.75000 

1.00000 

843 

4.90000 

1 .33333 

751 

4.32854 

1.04438 

844 

4.95000 

1 .33333 

752 

4.39569 

1.04438 

845 

5.00000 

1.33333 

753 

4.46283 

1.04438 

846 

5.05000 

1.33333 

754 

4.52998 

1.04438 

847 

5.10000 

1.33333 

755 

4.59713 

1.04438 

848 

5.15000 

1 .33333 

756 

4.66427 

1.04438 

849 

5.20000 

1 .33333 

757 

4.73142 

1.04438 

850 

5.25000 

1.33333 

758 

4.79856 

1.04438 

851 

5.30000 

1 .33333 

759 

4.86571 

1.0U38 

852 

5.35000 

1.33333 

760 

4.93285 

1.04438 

853 

5.40000 

1 .33333 

761 

5.00000 

1.04438 

854 

5.45000 

1.33333 

762 

5.06714 

1.04439 

855 

5.50000 

1 .33333 

763 

5.13429 

1.04439 

856 

4.50000 

1.41667 

764 

5.20144 

1.04439 

857 

4.55000 

1.41667 

765 

5.26858 

1.04439 

858 

4.60000 

1.41667 

766 

5.33573 

1.04439 

859 

4.65000 

1.41667 

767 

5.40287 

1.04439 

860 

4.70000 

1.41667 

768 

5.47002 

1.04439 

861 

4.75000 

1.41667 

769 

5.53716 

1.04439 

862 

4.80000 

1.41667 

770 

5.60431 

1.04439 

863 

4.85000 

1.41667 

771 

5.67145 

1.04439 

864 

4.90000 

1.41667 

772 

4.39783 

1.10216 

865 

4.95000 

1.41667 

773 

4.45805 

1.10216 

866 

5.00000 

1.41667 

774 

4.51827 

1.10216 

867 

5.05000 

1.41667 

775 

4.57848 

1.10216 

868 

5.10000 

1.41667 

776 

4.63870 

1.10216 

869 

5.15000 

1.41667 

777 

4.69892 

1.10216 

870 

5.20000 

1.41667 

778 

4.75913 

1.10216 

871 

5.25000 

1.41667 

779 

4.81935 

1.10216 

872 

5.30000 

1.41667 

780 

4.87956 

1.10216 

873 

5.35000 

1.41667 

781 

4.93978 

1.10216 

874 

5.40000 

1.41667 

782 

5.00000 

1.10216 

875 

5.45000 

1.41667 

783 

5.06021 

1.10216 

876 

5.50000 

1.41667 

784 

5.12043 

1.10216 

877 

4.50000 

1.50000 

785 

5.18065 

1.10216 

878 

4.55000 

1.50000 

786 

5.24086 

1.10216 

879 

4.60000 

1.50000 

787 

5.30108 

1.10216 

880 

4.65000 

1.50000 

788 

5.36130 

1.10216 

881 

4.70000 

1.50000 

789 

5.42151 

1.10216 

882 

4.75000 

1.50000 

790 

5.48173 

1.10216 

883 

4.80000 

1.50000 

791 

5.54195 

1.10216 

884 

4.85000 

1.50000 

792 

5.60216 

1.10216 

885 

4.90000 

1.50000 

793 

4.45561 

1.17145 

886 

4.95000 

1.50000 

794 

4.51005 

1.17145 

887 

5.00000 

1.50000 

795 

4.56449 

1.17145 

888 

5.05000 

1.50000 

796 

4.61893 

1.17145 

889 

5.10000 

1.50000 

797 

4.67337 

1.17145 

890 

5.15000 

1.50000 

798 

4.72781 

1.17145 

891 

5.20000 

1.50000 

799 

4.78225 

1.17145 

892 

5.25000 

1.50000 

800 

4.83668 

1.17145 

893 

5.30000 

1.50000 

801 

4.89112 

1.17145 

894 

5.35000 

1.50000 

802 

4.94556 

1.17145 

895 

5.40000 

1.50000 

803 

5.00000 

1.17145 

896 

5.45000 

1.50000 

804 

5.05444 

1.17145 

897 

5.50000 

1.50000 

805 

5.10888 

1.17145 

898 

5.75000 

0.500000 

806 

5.16331 

1.17146 

899 

5.86250 

0.500000 

807 

5.21775 

1.17146 

900 

5.97500 

0.500000 

808 

5.27219 

1.17146 

901 

6.08750 

0.500000 

809 

5.32663 

1.17146 

902 

6.20000 

0.500000 

810 

5.38107 

1.17146 

903 

6.31250 

0.500000 

811 

5.43551 

1.17146 

904 

6.42500 

0.500000 

812 

5.48995 

1.17146 

905 

6.53750 

0.500000 

813 

5.54438 

1.17146 

906 

6.65000 

0.500000 

814 

4.50000 

1.25000 

907 

6.76250 

0.500000 

815 

4.55000 

1.25000 

908 

6.87500 

0.500000 

816 

4.60000 

1.25000 

909 

6.98750 

0.500000 

817 

4.65000 

1.25000 

910 

7.10000 

0.500000 

818 

4.70000 

1.25000 

911 

7.21250 

0.500000 

819 

4.75000 

1.25000 

912 

7.32500 

0.500000 

820 

4.80000 

1.25000 

913 

7.43750 

0.500000 

821 

4.85000 

1 .25000 

914 

7.55000 

0.500000 

822 

4.90000 

1.25000 

915 

7.66250 

0.500000 

823 

4.95000 

1.25000 

916 

7.77500 

0.500000 

824 

5.00000 

1.25000 

917 

7.88750 

0.500000 

825 

5.05000 

1.25000 

918 

8.00000 

0.500000 

826 

5.10000 

1.25000 

919 

5.75000 

0.600000 

827 

5.15000 

1.25000 

920 

5.86250 

0.597000 

828 

5.20000 

1.25000 

921 

5.97500 

0.594000 

829 

5.25000 

1.25000 

922 

6.08750 

0.591000 

830 

5.30000 

1.25000 

923 

6.20000 

0.588000 

831 

5.35000 

1 .25000 

924 

6.31250 

0.585000 

832 

5.40000 

1.25000 

925 

6.42500 

0.582000 

833 

5.45000 

1.25000 

926 

6.53750 

0.579000 

834 

5.50000 

1.25000 

927 

6.65000 

0.576000 

835 

4.50000 

1.33333 

928 

6.76250 

0.573000 

836 

4.55000 

1 .33333 

929 

6.87500 

0.570000 

837 

4.60000 

1 .33333 

930 

6.98750 

0.567000 

20 


931 

7.10000 

0.564000 

1024 

8.00000 

932 

7.21250 

0.561000 

1025 

8.10364 

933 

7.32500 

0.558000 

1026 

8.20617 

934 

7.43750 

0.555000 

1027 

8.30647 

935 

7.5500) 

0.552000 

1028 

8.40348 

936 

7.66250 

0.549000 

1029 

8.50000 

937 

7.77500 

0.546000 

1030 

8.58333 

938 

7.88750 

0.543000 

1031 

8.66667 

939 

8.00000 

0.540000 

1032 

8.75000 

940 

5.75000 

0.700000 

1033 

8.00000 

941 

5.86250 

0.694000 

1034 

8.08581 

942 

5.97500 

0.688000 

1035 

8.17163 

943 

6.08750 

0.682000 

1036 

8.25744 

944 

6.20000 

0.676000 

1037 

8.34325 

945 

6.31250 

0.670000 

1038 

8.42906 

946 

6.42500 

0.664000 

1039 

8.51488 

947 

6.53750 

0.658000 

1040 

8.60069 

948 

6.65000 

0.652000 

1041 

8.68650 

949 

6.76250 

0.646000 

1042 

8.00000 

950 

6.87500 

0.640000 

1043 

8.07858 

951 

6.98750 

0.634000 

1044 

8.15715 

952 

7.10000 

0.628000 

1045 

8.23573 

953 

7.21250 

0.622000 

1046 

8.31431 

954 

7.32500 

0.616000 

1047 

8.39289 

955 

7.43750 

0.610000 

1048 

8.47146 

956 

7.55000 

0.604000 

1049 

8.55004 

957 

7.66250 

0.598000 

1050 

8.62862 

958 

7.77500 

0.592000 

1051 

8.00000 

959 

7.88750 

0.586000 

1052 

8.07220 

960 

8.00000 

0.580000 

1053 

8.14439 

961 

5.75000 

0.800000 

1054 

8.21658 

962 

5.86250 

0.791000 

1055 

8.28878 

963 

5.97500 

0.782000 

1056 

8.36097 

964 

6.08750 

0.773000 

1057 

8.43317 

965 

6.20000 

0.764000 

1058 

8.50536 

966 

6.31250 

0.755000 

1059 

8.57756 

967 

6.42500 

0.746000 

1060 

8.00000 

968 

6.53750 

0.737000 

1061 

8.06680 

969 

6.65000 

0.728000 

1062 

8.13360 

970 

6.76250 

0.719000 

1063 

8.20039 

971 

6.87500 

0.710000 

1064 

8.26719 

972 

6.98750 

0.701000 

1065 

8.33399 

973 

7.10000 

0.692000 

1066 

8.40079 

974 

7.21250 

0.683000 

1067 

8.46759 

975 

7.32500 

0.674000 

1068 

8.53438 

976 

7.43750 

0.665000 

1069 

8.00000 

977 

7.55000 

0.656000 

1070 

8.06250 

978 

7.66250 

0.647000 

1071 

8.12500 

979 

7.77500 

0.638000 

1072 

8.18750 

980 

7.88750 

0.629000 

1073 

8.25000 

981 

8.00000 

0.620000 

1074 

8.31250 

982 

5.75000 

0.900000 

1075 

8.37500 

983 

5.86250 

0.888000 

1076 

8.43750 

984 

5.97500 

0.876000 

1077 

8.50000 

985 

6.08750 

0.864000 

1078 

10.0000 

986 

6.20000 

0.852000 

1079 

10.0000 

987 

6.31250 

0.840000 

1080 

10.0000 

988 

6.42500 

0.828000 

1081 

10.0000 

989 

6.53750 

0.816000 

1082 

9.86667 

990 

6.65000 

0.804000 

1083 

9.85556 

991 

6.76250 

0.792000 

1084 

9.84444 

992 

6.87500 

0.780000 

1085 

9.83333 

993 

6.98750 

0.768000 

1086 

9.73333 

994 

7.10000 

0.756000 

1087 

9.71111 

995 

7.21250 

0.744000 

1088 

9.68889 

996 

7.32500 

0.732000 

1089 

9.66667 

997 

7.43750 

0.720000 

1090 

9.60000 

998 

7.55000 

0.708000 

1091 

9.56667 

999 

7.66250 

0.696000 

1092 

9.53333 

1000 

7.77500 

0.684000 

1093 

9.50000 

1001 

7.88750 

0.672000 

1094 

9.51500 

1002 

8.00000 

0.660000 

1095 

9.47667 

1003 

5.75000 

1.00000 

1096 

9.43833 

1004 

5.86250 

0.985000 

1097 

9.40000 

1005 

5.97500 

0.970000 

1098 

9.43000 

1006 

6.08750 

0.955000 

1099 

9.38667 

1007 

6.20<K)0 

0.940000 

1100 
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0.880000 

1104 

9.24833 

1012 

6.76250 
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0.745000 

1113 

9.00000 

1021 

7.77500 

0.730000 

1114 

9.09000 

1022 

7.88750 

0.715000 

1115 

9.02667 

1023 

8.00000 

0.700000 

1116 

8.96333 

0.700000 

1117 

8.90000 

1.09500 

0.705385 

1118 

9.00500 

0.930000 

0.721483 

1119 

8.93667 

0.962500 

0.748120 

1120 

8.86833 

0.995000 

0.785010 

1121 

8.80000 

1.02750 

0.825000 

1122 

8.92000 

0.870000 

0.800000 

1123 

8.84667 

0.900000 

0.775000 

1124 

8.77333 

0.930000 

0.750000 

1125 

8.70000 

0.960000 

0.660000 

1126 

8.83500 

0.810000 

0.664308 

1127 

8.75667 

0.837500 

0.677186 

1128 

8.67833 

0.865000 

0.698496 

1129 

8.60000 

0.892500 

0.728008 

1130 

8.75000 

0.750000 

0.760000 

1131 

8.66667 

0.775000 

0.740000 

1132 

8.58333 

0.800000 

0.720000 

1133 

8.50000 

0.825000 

0.71S616 
0.620000 
0.623231 
0.632890 
0.648872 
0.671006 
0.695000 
0.680000 
0.665000 
0.672444 
0.580000 
0.582154 
0.588593 
0.599248 
0.614004 
0.630000 
0.620000 
0.610000 
0.621383 
0.540000 
0.541077 
0.544297 
0.549624 
0.557002 
0.565000 
0.560000 
0.555000 
0.563500 
0.500000 
0.500000 
0.500000 
0.500000 
0.500000 
0.500000 
0.500000 
0.500000 
0.500000 

1.35000 

1.40000 

1.45000 

1.50000 

1.35000 

1.40000 

1.45000 

1.50000 

1.35000 

1.40000 

1.45000 

1.50000 

1.35000 

1.40000 

1.45000 

1.50000 
1.29000 
1.33750 
1.38500 
1.43250 
1.23000 
1.27500 
1.32000 
1.36500 
1.17000 
1.21250 
1.25500 
1.29750 
1.11000 
1.15000 
1.19000 
1 .23000 
1 .05000 
1.08750 
1.12500 
1.16250 
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1.06000 
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MESH  GENERATOR  ••  AUTOMATIC  NESM  GENERATION  FOR  FILE 
‘MESH. IN'  -•  MESHES  ALL  AREAS  IN  THE  FILE 

IN  THIS  CODE  'X'  COORDINATES  ARE  THE  *Z'  AXIS 
.  'Y*  COORDINATES  ARE  THE  'R'  AXIS 


CHARACTER*!  RIN 
DIMENSION  NE(8) 

WRITE  (»,*) 

WRITE  (*  *)  'CHOOSE:' 

WRITE  (*!*)  '1  •  MESH  NEW  'MCSH.IN*  FILE' 

WRITE  (*,*)  '2  -  MERGE  PERIMETER  NODES' 

WRITE  (*.*)  '3  -  CONDENSE  OUT  UNUSED  NODE  NUMBERS' 
WRITE  (*,*)  '*  -  (NOT  USED)' 

WRITE  (*.*)  '5  •  END  (WRITE  OUTPUT  FILES)' 

WRITE  (*,*)  '  •  INCLUDING  F.E.  COMPATIBLE  INPUT  FILE* 

WRITE  (* ,*) 

READ  (*  *)  ICHOICE 
WRITE  (*  *) 

IF  (ICHoicE.EQ.2)  GO  TO  1500 
IF  (ICHOICE. EQ.3)  GO  TO  2500 
IF  (ICHOICE. EQ. 4)  GO  TO  3000 
IF  (ICHOICE. EQ. 5)  GO  TO  4000 

OPEN  (1,FILE«'MESH.IN>,STATUS«'K}L0".  FORN=»FORMATTEO», 
RECL«1 ) 


READ  ENTIRE  INPUT  FILE  AND  CREATE  A  FORMATTED  SCRATCH  FILE 
POINTS 

OPEN  (4,FILE«'INPUT'. ACCESS-'V I RECT* . FORM^" FORMATTED" , 
RECL-80 . STATUS»"NEU" ) 

OPEN  ( 7, F I LE» ' PLOT .IN', STATUS«"NEW" ) 

NUMP>0 

FORMAT  (1X,I12,F15.6,F15.6) 

FORMAT  (1X. 5(112), 15) 

NUML>0 

NUMS«0 

NUMA«0 

FIND  ALL  POINTS 

FORMAT  (A1,I5.F15.6,F1S.6) 

READ  (1,*)  RIN 
IF  (RIN.EQ.'*')  GO  TO  30 
IF  (RIN.NE.'P')  GO  TO  10 
WRITE  (*,*)  RIN 
BACKSPACE  (1) 

X>0.0 

Y*0,0 

N^O 


READ  (1,*)  RIN.N,X,Y 

WRITE  (*,*)  RIN,N,X,Y 

IF  (N.GT.NUMP)  NUMP«N 

WRITE  (4,20,RECsN)  RIN,N,X,Y 

GO  TO  10 

FIND  ALL  LINES 

REWIND  (1) 

FORMAT  (A1,S(I5)) 

FORMAT  (A1,7(I5)) 

READ  (1,*)  RIN 
IF  (RIN.EQ.'*')  GO  TO  50 
IF  (RIN.NE.'L')  GO  TO  45 
WRITE  (*,*)  RIN 
BACKSPACE  (1) 

N>0 

IC>0 

READ  (1,*)  RIN,N,I1.I2.N0,IC 
WRITE  (*,*)  RIN,N,I1,I2,N0.IC 
IF  (N.GT.NUML)  NUML>N 
IREC>NUMP'»N 

WRITE  (4,40,RECsIREC)  RIN,N.I1,I2,ND.IC 
GO  TO  45 

FIND  ALL  SPLINES 
REWIND  (1) 

FORMAT  (A1,6(I5)) 

READ  (1,*)  RIN 
IF  (RIN.EQ.'*')  GO  TO  70 
IF  (RIN.NE.'S')  GO  TO  65 
WRITE  (* ,*)  RIN 
BACKSPACE  (1) 

N«0 

L1»0 


L2«0 

L3-0 

L4-0 

L5-0 

READ  (1,*)  RIN,N,L1,L2.L3.L4,L5 
WRITE  (^.*)  RIN,N,L1,L2,lS.L4,L5 
IF  (N.GT.NUMS)  NUMS«N 
IREC>NUMP«NUML-HI 

WRITE  (4,60,REC>IREC)  RIN,N,L1,L2.l3.L4,LS 

GO  TO  65 

READ  ALL  AREAS 

REWIND  (1) 

READ  (1,*)  RIN 


33 


c 

c 

151 

150 


155 


160 

C 

C 

C 

c 

C165 

C 

C 

C 

C 

Cl  70 
C 


c 


c 


c 

c 


c 

c 

c 

c 

c 

c 


200 

250 

c 

c 

c 


IF  (RIN.Efl.'*')  GO  TO  ISO 
IF  (RIN.ME.'A*)  GO  TO  7S 
WRITE  (*,*)  RIN 
BACKSPACE  (1) 

N«0 

J1-0 

J2«0 

J3-0 

J4»0 

READ  (1,*)  RIN,N,J1,J2.J3.J4,NTYPE,NMATER1AL 
WRITE  (*.*)  RIN,N,Jl,J2,J3,j4,NTyPE.NMATERIAL 
IF  (N.GT.NUMA)  NUMA«N 
I  REC<NUMP«IHJML>NUMS'»N 

WRITE  (4,41,REC«IREC)  RIH.M, J1. J2, J3, J4,MTYPE,NNATERIAL 
GO  TO  75 


WRITE  A  POINT  AND  LINE  PLOT  INPUT  FILE  FOR  MESHPLOT 
FORMAT  (A1.I12,F15.6.F15.6.I9.I9) 

WRITE  (7,*)  ‘N'.NUMP.HUML.O.O.O 

DO  1SS  I«1,NUMP 

IREC>I 

READ  (4.20.REC>IREC)  RIN.N.X.Y 

WRITE  (7,1S1)  RIN.N.X.Y. 6,0 

CONTINUE 

DO  160  I>1,NUML 

IRECaNUMP^I 

READ  (4>0,REC«IREC)  RIN. N. II,  12. 13, 14 

WRITE  (7,*)  RIN. N. 11, 12,13, 14 

CONTINUE 

DO  165  1-1, HUMS 

IREC>NUMP>NUML«I 

READ  (4.60,REC«IREC)  RIN,N,L1,L2,L3,L4,L5 
WRITE  (7.*)  RIN,N.L1.L2,L3.L4,L5 
CONTINUE 
00  170  l«1,NUNA 
I  REC>NUMP>NUNL'*'NUMS*  I 

READ  (4.41,REC«IREC}  RIN.N.JI, J2, J3, J4, IM.NM 
WRITE  (7  •}  RIN,N,J1,J2.J3,J4 
CONTINUE 
CLOSE  <7) 

STARTING  NODE  NUMBER 
IN>0 

CLOSE  (1) 

OPEN  ( 1 ,  F I  LE« '  MERGE  • , ACCESS«“OIRECT'',  FORM="FORMATTED" , 
RECLaS.STATUS-xNEW") 

OPEN  (2,FILE»'NOOESL*,ACCESS»"DIRECT'',FORM="FORMATTED», 
RECL*43.STATUS«'‘NEW'‘) 

STARTING  ELEMENT  HUMBER 
IE«0 

OPEN  (3,FILE»'ELEML',ACCESS»"OIRECT'',FORM="FORMATTED", 
RECL«66,STATUS>"NEU») 


00  1000  IA«1,I 
I  REC«NUMP'»NUML'»NUMS'»  I A 
N«0 
JPl-0 
JP2-0 
JP3»0 
JP4-0 
MTYPE=0 

READ  (4,41,REC-IREC)  RIN,N,JP1,JP2, JP3, JP4,HTYPE,NMAT 
IF  (N.EQ.O)  GO  TO  1000 

WRITE  (*,*)  RIN,M,JP1,JP2,JP3,JP4,HTYPE,MMAT 

N01«0 

N02«0 


WRITE  <*,*)  'CALL  EOGEMESH  K»1' 

FINO  ANO  MESH  LINE.  RADIUS,  OR  SPLINE  CONNECTING  JP1  TO  JP2 
CALL  EOGEMESH  <1, JP1, JP2,N01,N02.IN,NUMP,NUML,NUMS> 

WRITE  (*,*)  'CALL  EOGEMESH  K>2' 

FIM)  ANO  MESH  LINE,  RADIUS,  OR  SPLINE  CONNECTING  JP2  TO  JP3 
CALL  EOGEMESN  (2, JP2,JP3,N01,ND2, IN.NUMP.NUML.NUMS) 

FINO  ANO  MESH  LINE.  RADIUS,  OR  SPLINE  CONNECTING  JP3  TO  JP4 
CALL  EOGEMESN  (3, JP4,JP3.ND1,ND2. IN.NUMP.NUML.NUMS) 

FINO  ANO  MESH  LINE.  RADIUS,  6r  SPLINE  CONNECTING  JP4  TO  JP1 
CALL  EOGEMESH  <4, JP1,JP4,H01,ND2. IN.NUMP.NUML.NUMS) 

ASSIGN  NODES  TO  ELEMENTS  FOR  THIS  AREA 

00  250  KE>1,N02 
DO  200  JEal.NOI 
lE-IE-'^l 

H1«IN+JE+((N01*1)*(KE-1)) 

N2-HU1 

N3aN2'»N0U1 

N4>N3-1 

WRITE  (3.8,REC>IE)  1E,H1,N2,N3,N4,NMAT 

CONTINUE 

CONTINUE 


FINO  X.Y  COORDINATES  FOR  INTERIOR  NODES  FOR  THIS  AREA 
NNalN^UNOUl 


34 


300 

400 

1000 

1001 


C 

C 

1500 


c 


1510 

C 

C 

C 

1610 

1611 

1615 


C 

1620 


C 

C 


00  400  l>1,(N02-1) 

NEI-IN^IHiDUdHOI)*! 

NS1>NE1*N01 

READ  (2,7,REC«NE1:  NX.XEI.YEI 
WRITE  (*  *)  NX.XEI.YEI 
READ  (2.7,REC>NS1}  NX.XSI.YSI 
0X«(XE1-XS1}/NO1 
00  300  J«1.(N01-1} 

IF  (NTYPE.E0.2)  THEN 

NE2>IN«(N0U1)*(N02«1}-N0UJ 

NS2>IN4'UJ 

READ  (2.7.REC-NE2)  NY.XE2.YE2 

READ  (2.7.REC-NS2)  NY.XS2.YS2 

0Y>(YE2-YS2)/N02 

Y«YS2>0Y*I 

ENOIF 

IF  (NTYPE.EQ.I)  THEN 

0Y«(YE1-YS1)/N01 

Y»YS1*OY*J 

ENOIF 

NN-NN+1 

XsXS1'H>X*J 

WRITE  (2.7.REC«NN)  NN.X.Y 

CONTINUE 

NN«NN>2 

CONTINUE 

IN-IN»(N01«1)*(N02«1) 

CONTINUE 
FORMAT  (15) 

WRITE  (1.1001)  0 
CLOSE  (1) 

'CLOSE  (2) 

CLOSE  (3) 

CLOSE  (4) 

CLOSE  (5) 

CLOSE  (6) 

OPEN  ( 1 . F I LE- *  TOTALS ‘ . STATUS>"NEW» ) 

WRITE  (1.*)  IN.IE 
CLOSE  (1) 

GO  TO  1 

MERGE  NODES  ROUTINE 
OPEN  d.FILE-'TOTALS') 

READ  (1  *)  IN.IE 
CLOSE  (1) 

OPEN  (2.FILE»'NOO£SL'.ACCESS»“OIRECT".FORM=''FORMATTED", 
RECL«43) 

OPEN  (3.FILE><ELEML'. ACCESS»"0 1 RECT . FORM»''FORMATTEO" , 
RECL>66) 

OPEN  (4 .  F I  LE« '  MERGE ' .  ACCESS-'*0 1  RECT" .  FORM*"  FORMATTED" , 
RECL*5) 

CREATE  A  NODE-ELEMENT  FILE 

OPEN  (5. FILE* 'NOEL*. ACCESS>«0 I RECT" . FORM*" FORMATTED" , 
RECL*45) 

FORMAT  (9(15)) 

CALL  NOEL  (IN.IE) 

MERGE  OUT  NODES  ON  COMMON  PERIMETER  LINES 

FORMAT  (A2) 

FORMAT  (15) 

JREC-1 

KREC>JREC 

READ  (4.1611.REC>KREC)  ISTAR 
WRITE  (*  •)  'ISTAR*. ISTAR 
IF  (ISTAR.EQ.O)  GO  TO  2400 
READ  (4,1611.REC-KREC)  LI 
WRITE  (*.*)  'LI*'. Li 
KREC*KREC'*'1 

READ  (4.1611.REC*KREC)  NNL1 

KREC*KREC«NNLU1 

JREC*KREC 

IF  (L1.EQ.-1)  GO  TO  1615 

FIND  IF  LINE  Li  REPEATS  IN  THIS  FILE  AND  MERGE  OUT  NODES 

READ  (4,1611.REC*KREC)  ISTAR 

WHITE  (*.•)  'ISTAR*. ISTAR 

IF  (ISTAR.EQ.O)  GO  TO  1615 

READ  (4.1611. REC*KREC)  L2 

WRITE  (*.*)  'L2*'.L2 

MREC*KREC 

KREC*KREC«^1 

READ  (4.1611,REC*KREC)  NNL2 
IF  (L1.NE.L2)  GO  TO  1650 
WE  HAVE  A  MATCH 
KREC*KREC«1 

WRITE  (4.1611.REC*MREC)  -1 

REPLACE  L2  NODES  WITH  Li  NOOES  IN  ELEMENT-NODE  FILE  'ELEML' 

LREC1*JREC-NNL1 

DO  1625  I*1,NNL1 

READ  (4.1611,REC*LREC1)  N1 

READ  (4,1611.REC*KREC)  N2 

LREC1*LRECU1 

KREC*KREC«1 


URITE  (*  *)  Ml  M2 

READ  (Sj510,Rk«N2)  JN2,NE<1>.ME(2).ME(3),ME(4),ME(S),ME(6), 

*  Me<7).N6<8) 

DO  1^4  J>1,8 

IF  (NE(J).EQ.O)  GO  TO  1624 
NREC>ME(J) 

REM)  (3.8,REC»«IREC}  N. J1. J2, J3, J4.N)UT 
IF  (M2.EQ.J1)  URITE  (3.8.REC-MREC)  N.MI. J2. J3. J4.MMAT 

IF  (N2.EQ.J2)  WRITE  (3,8,REC»MREC)  M,J1,N1,J3,J4,MMAT 

IF  (N2.EQ.J3)  URITE  (3,8,REC>NREC)  N.J1.J2,N1,J4,NMAT 

IF  (N2.EQ.J4)  WRITE  (3.8.REC-NREC)  M.JI. J2, J3,N1,NMAT 

C  REPUCE  M2  WITH  0  IM  'MOOESL'  FILE 

URITE  (2.7,REC>M2)  0.0,0 

1624  COMTIMJE 

1625  COMTIMUE 
GO  TO  1620 

1650  KREC«KREC«MML2«1 
1660  GO  TO  1620 
2400  COMTIMUE 
C  FIMISHEO 

CLOSE  (2) 

CLOSE  (3) 

CLOSE  (4) 

CLOSE  (5) 

GO  TO  1 
C 

C  CONOEMSE  MODE  MUH8ERS  ROUT (ME 

2500  OPEM  (1,F1LE«‘ TOTALS') 

READ  (1  *)  IM.IE 

OPEM  <2lFILE«*MOOESL',ACCESS>i"OIRECT",FORM="FORMATTEO", 

*  RECL>43) 

OPEM  <3,FILE»'ELEML‘.ACCESS*"OIRECT“.FORM="FORHATTED“, 

*  RECL>66) 

OPEM  <5,FIL£»*M0EL‘, ACCESS»“0 I RECT" , FORM*" FORMAT  TED " , 

*  RECL-45) 

CALL  MOEL  (IM.IE) 

ID>0 

DO  2600  IR*1,IM 
READ  (2,7,REC-IR)  M.X.Y 
URITE  (*,*)  M.X.Y 
IF  (M.EQ.O)  ID-IO^I 
IF  (M.EQ.O)  GO  TO  2600 
IF  (lO.EQ.O)  GO  TO  2600 
IU*IR*I0 

READ  (5,1510,REC«IR)  4M2,ME(1).ME(2),ME(3),ME(4),ME(5),ME(6), 

*  ME(7).NE(8) 

N>IR'I0 

URITE  (2,7,REC>IU)  M.X.Y 
URITE  (2,7,REC>1R)  0,0.0 
DO  2550  J«{,8 
IF  (ME(J).EO.O)  GO  TO  2550 
MREC*ME(J) 

READ  (3,8,REC>MREC)  JE, J1 . J2. J3, J4, JNAT 
IF  (IR.EQ.J1)  URITE  (3,8,REC«MREC)  JE.M, J2, J3, J4, JMAT 
IF  (IR.EQ.J2)  URITE  (3,8,REC-NREC)  JE, Jl,N. J3, J4, JMAT 
IF  (IR.EQ.J3)  URITE  (3,8,REC*MREC)  JE. Jl. J2,N, J4, JMAT 
IF  (IR.EQ.J4)  URITE  (3,8.REC>MREC)  JE, Jl. J2, J3,N, JMAT 
2550  COMTIMUE 
2600  COMTIMUE 
IM*IM-I0 
REUIMO  (1) 

URITE  (1,*)  IM.IE 
CLOSE  (1) 

CLOSE  (2) 

CLOSE  (3) 

CLOSE  (5) 

GO  TO  1 
C 

C  REMUMBER  MOOES,  AND  ELEMEMTS  (BANDUIDTH)  ROUTINE 

3000  COMTIMUE 

GO  TO  1 
C 

C  URITE  OUTPUTFILES  COMPATIBLE  WITH  HESHPLOT 

4000  OPEM  (5,FIlE»'MOOES',STATUS«"MEW") 

OPEM  (6, F I LE> 'ELEMEMTS • , STATUS>"NEU" ) 

OPEM  (2.  F I  LE> '  MOOESL  • , ACCESS»''OIRECT",  FORM*'' FORMATTED" , 

*  RECL>43) 

OPEM  (3,FILE*'ELEML' ,ACCESS»''OIRECT“,FORM*"FORMATTED", 

*  RECL«66) 

OPEM  (4,FILE-'FE.IM') 

4001  FORMAT  ('  ET, 1, 42,0,0,1 • ) 

4002  FORMAT  (•  EX,',I3 ,',') 

4003  FORMAT  ('  MUXY, ' , 13, ' , ' ) 

4004  FORMAT  ('  OEMS, ' , 13, ' , ' ) 

4005  FORMAT  ('  MAT, '.13) 

OPEM  (1.F1LE«'T&TALS') 

READ  (1,*)  IM.IE 
CLOSE  (1) 

DO  4010  1*1, IM 
READ  (2,7,REC«I)  M.X.Y 
WRITE  (5,*)  '  '.M.X.Y 
C  F.E.  USES  Y-Z  AMO  X*R 

URITE  (4,*)  '  M.'.M.'.'.Y.'.'.X 


4010  C0NTINU6 

MITE  (4,4001) 

HNMAT>0 

DO  4020  I>1,IE 

READ  (3,8,REC«I)  M.MI .N2,N3.N4,NMAT 
MITE  (6,^)  ■  <,M,N1,N2.N3,M4,tlNAT 
IF  (NHAT.EQ.NNMAT)  GO  TO  4015 
NNMAT-NNAT 

MITE  (4,4002)  NNMAT 
MITE  (4,4003)  NNMAT 
MITE  (4,4004)  NNMAT 
MITE  (4,4005)  NNMAT 

4015  MITE  (4,*)  •  E,',N1,',*,N2,',*,N3,',',N4 

4020  CONTINUE 

5000  ENO 

SUBROUTINE  NOEL  (IN. IE) 

1  FORMAT  (9(15)) 

2  FORMAT  (IX, 5(112), 15) 

00  10  1-1, IN 

10  MITE  (5,1,REC-I)  1 ,0,0,0,0,0,0,0.0 

DO  20  1-1, IE 

READ  (3,2,REC-I)  N.N1.N2.N3.N4.NMAT 
CALL  ENTER  (N1,I) 

CALL  ENTER  (N2,I) 

CALL  ENTER  (N3,I) 

CALL  ENTER  (N4,I) 

20  CONTINUE 

RETURN 
ENO 

SUBROUTINE  ENTER  (K.J) 

1  FORMAT  (9(15)) 

READ  (5,1.REC«K)  I .M1.M2,M3,M4,N5,M6,M7,M8 
IF  (M1.EO.0)  MITE  (5.1,REC-K)  I , J,N2,M3,M4,M5,M6.M7,H8 
IF  (M1.EO.O)  GO  TO  10 

IF  (M2.EO.O)  MITE  (5,1.REC>K)  I.M1.J,H3.M4,M5,M6,M7,MS 
IF  (M2.EQ.0)  GO  TO  10 

IF  (M3.EO.0)  MITE  (5,1,REC-K)  I,M1.H2.J.M4,H5,H6,N7.M8 
IF  (M3.EQ.0)  GO  TO  10 

IF  (M4.EQ.0)  MITE  (5,1,REC-X)  I,M1,M2,M3, J,HS,N6,H7.M8 
IF  (MA.EQ.O)  GO  TO  10 

IF  (M5.Ea.O)  MITE  (5,1,REC«K)  I,M1,M2,N3,H4,J,H6,H7,M8 
IF  (M5.EQ.0)  GO  TO  10 

IF  (M6.EQ.0)  MITE  (5,1,REC-K)  I,M1,M2,M3,M4,H5,J.M7,M8 
IF  (M6.EO.O}  GO  TO  10 

IF  (M7.EQ.0)  MITE  (5,1,REC>K)  I,M1,M2,N3,M4,M5,H6, J,M8 
IF  (M7.EO.O)  GO  TO  10 

IF  (M8.EQ.0)  MITE  (5.1,REC-K)  I,M1,M2,N3,M4,M5,H6,H7, J 
10  RETURN 

ENO 

SUBROUTINE  EOGEMESH  (K, JJ1 , JJ2,N01 ,N02, I IN,NUMP,NUML,NUHS} 

C  'K*  «  EDGE  NUMBER;  MIN'  -  ENDING  NODE  NUMBER  OF  LAST  AREA 

ruADArTFP-1  Pin 

DIMENSION  L(5),IP1(5),IP2(5),IND(5),IC(5) 

1  FORMAT  (15) 

2  FORMAT  (A1,I5.F15.6,F15.6) 

3  FORMAT  (A1) 

4  FORMAT  (A1,5(I5)) 

6  FORMAT  (A1,6(I5)) 

GO  TO  20 

C 

C  LOOK  FOR  A  SPLINE  CONNECTING  JJ1  AND  JJ2 

C 

7  IF  (NUMS.EQ.O)  GO  TO  100 
IREC-NUMP^MJML 

10  IREC-IRE01 

IF  (IREC.GT.(NUMI»NUML+NUMS))  GO  TO  100 
ICHECK-O 

READ  (4,6,REC-IREC)  RIN,N.L(1),L(2),L(3),L(4),L(5) 

MITE  (*,*)  RIN,N,L(1),L(2).L(3),L(4),L(5) 

JF-5 

IF  (L(2).Ea.O)  GO  TO  10 

IF  (L(5).EQ.O)  JF-4 

IF  (L(4).Ea.O)  JF-3 

IF  (L(3).E0.0)  JF-2 

16  JREC-NUMP«L(1) 

READ  (4,4,REC-JREC)  RIN,L(1). IP1(1), IP2(1 ). IND(1 ), IC(1 ) 

MITE  (*,*)  RIN,L(1),IP1(1),IP2(1),IND(1),IC(1) 

IF  (IP2(1).EQ.JJ1)  THEN 

IRP-IPKD 

IP1(1)»!P2(1) 

IP2(1)-IRP 

ENOIF 

IF  (IPI(I).EO.JJI)  GO  TO  17  37 

IF  (ICHECK.EQ.I)  GO  TO  10 
C  REVERSE  SPLINE  ORDER  AND  THEN  CHECK 

L1«L(5) 

L2-L(4) 

L3-L(3) 


rtoo 


L4>L(2) 

L5»L{1) 

IF  (JF.EO.S)  THEM 

L(1)«L1 

L(2)»L2 

L(3)»L3 

L(4)>L4 

L(5)-t5 

ENOIF 

IF  (JF.Ea.4}  THEN 

L(1)»L2 

L<2)«L3 

L(3)»L4 

L(4)>LS 

L(5)«L1 

ENOIF 

IF  (JF.EO.S)  THEN 

L(1)»L3 

L(2)>L4 

L(3)»l.5 

L(4)«L1 

L(5)«L2 

ENOIF 

IF  (JF.EQ.2)  THEN 

L(1)«L4 

L(2)«L5 

L(3)«L1 

L(4)«L2 

L<5)»L3 

ENOIF 

I CHECK- 1 

(»  TO  16 

17  JREC-NUMP«L(JF) 

REAO  (4,4,REC>JREC)  RIN.L(JF). IPI(JF). 1P2(JF), IND(JF), tC(JF) 
WRITE  (*.*)  RIN.L(JF},IP1(JF).IP2(JF).IN0(JF),1C(JF) 

IF  (IP1(JF).EQ.JJ2)  THEN 
IRP>IP2(JF) 

IP2(JF)-IP1(JF) 

IP1(JF)>IRP 

ENOIF 

IF  (IP2(JF).NE.JJ2)  (U  TO  10 
WRITE  (*,*)  'ITS  A  SPLINE' 

C  IT'S  A  SPLINE  ***• 

C  REAO  THE  INTERIOR  LINES 

IF  (JF.EQ.2)  GO  TO  13 
00  12  I«1,JF 
IF  (I.E0.1)  GO  TO  11 
IF  (I.EQ.JF)  GO  TO  11 
JREC>NUNP>L(I} 

REAO  (4,4.REC-JREC)  RIN,L(I),IP1(I>,1P2(I),IN0(I),IC(I) 
WRITE  (*  *)  RIN,L(I),IPUI),IP2(I),IND(I),IC(I) 

IF  (IP1(I).NE.IP2(I*b)  then 
IRP«IP1(I) 

IP1(I)>IP2(I) 

IP2(I)-IRP 

ENOIF 

r  CONTINUE 

12  CONTINUE 

13  CONTINUE 
NO-0 

DO  14  1-1, JF 

14  NO-NO«INO(I) 

IF  (K.EQ.D  THEN 
N01-M) 

NS-IIN^I 

NIN-1 

ENOIF 

IF  (K.Ea.2)  THEN 
N02-N0 

HS-IIN-'^(1'»N01) 

NIN-UN01 

ENOIF 

IF  (K.EQ.3)  THEN 

NS-IIN4'(1>N01)*N02-'’1 

NIN-1 

ENOIF 

IF  (K.Ea.4)  THEN 
NS-IIN'''1 
NIN-UN01 
ENOIF 

00  15  1-1, JF 

WRITE  (*,*)  'CALL  MESH  (SPLINE  LINE)' 

WRITE  (1,1)  L(() 

WRITE  (1,1)  IN0(I)>1 

CALL  MESH  (NS,NIH, IPKI ), IP2(I ), INDd ), ICd )) 
NS-NS-fNIN^INOd) 

15  CONTINUE 
GO  TO  100 
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LOOK  FOR  A  LINE  OR  RADIUS 

IREC-NUW> 

IREC-IREC+1 


••  ■ 


% 


IF  (IREC.GT.dRJMP^HUNL))  GO  TO  7 

READ  (i,4.REC«IREC)  RIH,N.IP1(t),lP2(1),INO(1).lC(1) 

WRITE  (%*)  RlN,N,IP1(1),fP2<1),INO(1),IC(1) 

IF  (IPKb.EG.JJZ)  THEM 
IRP-IPKI) 

IPT(1)-IP2(T> 

IP2(1)«IRP 

ENOIF 

IF  (IPKD.NE.JJI)  GO  TO  2S 
IF  (IP2(1).NE.JJ2)  GO  TO  2S 
IF  (K.EQ.1)  THEN 
N01«IN0(1) 

HS>1IH-»I 

NIH-1 

ENOIF 

IF  (K.Ea.2)  THEN 
N02>IN0(1) 

NS>IIN«UN01 

NIN-1+N01 

ENOIF 

IF  (K.EQ.3)  THEN 

NS>11N«<UN01)*N02+1 

NIN«1 

ENOIF 

IF  (K.E0.4)  THEN 
NS>IIN«1 
NIN-UN01 
ENOIF 

WRITE  (1.1)  N 

WRITE  (1.1)  INOlD'fl 

WRITE  (*.*)  'CALL  HESH  LINE  OR  RADIUS' 

CALL  NESH  (NS.NIN. IP1(1).IP2(1). INO(1 ), IC(1 )) 

100  RETURN 
ENO 
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SUBROUTINE  MESH  (NS.NIN. II. I2.N0.JC) 

CHARACTER*!  RIN 

FORMAT  (A1.IS.F15.6.F1S.6) 

FORMAT  (1X.I12.F15.6.F15.6) 

FORMAT  (IS) 

IREC-I1 

READ  (4,1.REC>IREC)  RIN.N.XI.YI 
WRITE  (*.*)  RIN.N.XI.YI 
IREC-I2 

READ  (4.1.REC*IREC)  RIN.N.X2.Y2 

IF  (JC.NE.O)  THEN 

IREC«JC 

READ  (4.1.REC>IREC)  RIN.N.XC.YC 

ENOIF 

I«NS 

IF  (JC.NE.O)  GO  TO  20 

0Xa(X2-X1)/N0 

XS-X1 

0Y»(Y2-Y1)/N0 

YS*Y1 

DO  10  JaO.NO 

XI-XS*OX*J 

YIaYS*OY*J 

WRITE  (2.2.REC-I)  I.XI.YI 

WRITE  (1.3)  I 

I«I*NIN 

CONTINUE 

00  TO  100 

R«((X1-XC)**2*(Y1-YC)**2)**.5 

0X1>A8S(X1*XC) 

0X2>ABS(X2-XC) 

DY1aABS(Y1*YC) 

0Y2«ABS(Y2'YC) 

PI»3.14159 

IF  (X1.GT.XC.AN0.Y1.GE.YC)  T1*ATAM(DY1/0X1) 

IF  (X2.GT.XC.AN0.Y2.GE.TC)  T2sATAN(0Y2/DX2) 

IF  (X1.GT.XC.AN0.Y1.LT.YC)  T1»2.0*PI-ATAN(0Y1/DX1) 

IF  (X2.GT.XC.AN0.Y2.LT.YC)  T2«2.0*PI'ATAN(DY2/DX2) 

IF  (X1.LT.XC.AN0.Y1.GE.TC)  T1«PI-ATAN(DY1/DX1) 

IF  (X2.LT.XC.AN0.Y2.GE.TC)  T2sPI-ATAN(DY2/DX2) 

IF  (X1.LT.XC.AN0.Y1.LT.YC)  TlaPI'fATANIOYI/DXI) 

IF  (X2.LT.XC.AN0.Y2.LT.YC)  T2«PI«ATAN(DY2/DX2) 

IF  (XI.EO.XC.ANO.YI.GT.YC)  T1«PI/2.0 
IF  (X2.E0.XC.AND.Y2.GT.YC)  T2aPI/2.0 
IF  (XI.EO.XC.ANO.YI.LT.YC)  TlaPI'^PI/E.O 
IF  (X2.EQ.XC.AN0.Y2.LT.YC)  T2-Pl«PI/2.0 
IF  (T1.GE.0.0.AN0.T1.LE. (PI/2.0). AND. T2.GE.(PI*PI/2.0}) 
T2«2.0*PI-T2 

IF  (T2.GE.0.0.AN0.T2.LE. (PI/2.0). AND. T1.GE.(PI'''PI/2.0)) 

T2«T2*2.0*PI 

DT«(T2-T1)/N0 

00  30  J>0,N0 

XIaXC*R*C&(T1*J*0T) 

YI»YC*R*SIN(T1+J*0T) 

WRITE  (2.2.REC-I)  I.XI.YI 

WRITE  (1.3)  I 

Ia)*NIN 

CONTINUE 
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RETURN 

END 


